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O , ABSTRACT 

' We calculate axisymmetric toroidal modes of magnetized neutron stars with a solid 

\^ , crust. We assume the interior of the star is threaded by a poloidal magnetic field that 

^SJ ' is continuous at the surface with the outside dipole field whose strength Bp at the 

magnetic pole is Bp ~ lO^^G. Since separation of variables is not possible for oscil- 
lations of magnetized stars, we employ finite series expansions of the perturbations 
Q^, using spherical harmonic functions to represent the angular dependence of the oscil- 

lation modes. For Bp ^ lO^^G, we find distinct mode sequences, in each of which the 
oscillation frequency of the toroidal mode slowly increases as the number of radial 
^ ' nodes of the eigenfunction increases. The frequency spectrum of the toroidal modes 

Ci . for Bp ~ lO^^G is largely different from that of the crustal toroidal modes of the non- 

magnetized model, although the frequency ranges are overlapped each other. This 
suggests that an interpretation of the observed QPOs based on the magnetic toroidal 
^ ' modes may be possible if the field strength of the star is as strong as Bp ~ lO^^G. 

\o ' 

QQ I Key words: stars: neutron - stars: oscillations - stars : magnetic fields 

d 



o 



o 



1 INTRODUCTION 



Recent discovery of quasi-periodic oscillations (QPOs) of magnetar candidates is one of the observational manifestations of 
global oscillations of neutron stars. Israel et al (2005) detected QPOs of frequencies ~ 18, ~ 30 and ~92.5Hz in the tail of 
the SGR 1806-20 hyperflare observed December 2004, and suggested that the 30Hz and 92.5Hz QPOs could be caused by 
k>( ^ seismic vibrations of the neutron star crust (see, e.g., Duncan 1998). Later on, in the hyperflare of SGR 1900-1-14 detected 
' August 1998, Strohmayer & Watts (2005) found QPOs of frequencies 28, 53.5, 84, and 155 Hz, and claimed that the QPOs 
could be identified with the low / fundamental toroidal torsional modes of the solid crust of the neutron star. These recent 
discoveries of QPOs in the giant flares of Soft Gamma-Ray Repeaters SGR 1806-20 (Israel et al 2005, Watts & Strohmayer 
2006, Strohmayer & Watts 2006) and SGR 1900-1-14 (Strohmayer & Watts 2005) have made promising asteroseismology for 
magnetars, neutron stars with an extremely strong magnetic field (see, e.g., Woods &i Thompson 2006 for a review of SGRs). 

It is currently common to identify these QPOs with seismic vibrations caused by crustal toroidal modes of the neutron 
stars, since the frequency range of the modes overlap that of the observed QPOs and from the energetics point of view the 
crustal toroidal modes would be most easily excited to observable amplitudes by spending a least amount of available energies, 
for example, those released in magnetic field restructuring (e.g., Duncan 1998). Although the interpretation based on crustal 
torsional modes looks promising, we need detailed theoretical analyses of oscillations of magnetized neutron stars so that 
we could get information of physical conditions of the stars through the confrontation between theoretical modelings and 
observations. This is particularly true for high frequency QPOs (e.g., 625Hz QPO, Watts & Strohmayer 2006; 1835Hz QPO 
and less significant QPOs at 720 and 2384 Hz in SGR 1806-20, Strohmayer & Watts 2006), since there exist classes of modes 
other than the crustal toroidal modes that can generate the frequencies observed. 

The presence of a magnetic field makes it possible for toroidal modes to exist in a fluid star even without rotation as 
does the presence of the shear modulus in the solid crust. In this paper we are interested in axisymmetric toroidal modes 
since axisymmetirc toroidal and spheroidal modes are decoupled for a poloidal held when the star is non-rotating. For non- 
axisymmetric modes, the toroidal and spheroidal components are coupled even without rotation and hence the modal analyses 
of magnetized stars would be much more complicated. 
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Theoretical calculations of toroidal modes of strongly magnetized neutron stars have been carried out by several authors, 
including Piro (2005), Glampedakis et al (2006), Sotani et al (2006, 2007), and Lee (2007). The analyses by Piro (2005), 
Glampedakis et al (2006), and Lcc (2007) assume Newtonian gravity, while those by Sotani ct al (2006, 2007) use general 
relativistic formulation. Although the studies by Piro (2006) and Lee (2007) ignore the effects of magnetic fields in the fluid 
core, those by Glampedakis et al (2006) and Sotani et al (2006, 2007) consider magnetic waves propagating in the fluid core, 
assuming the core is threaded by a magnetic field of substantial strength. Besides the differences mentioned above, most of 
the authors except Lee (2007) represent the angular dependence of the oscillations by a single spherical harmonic function 
Yi^{9,(j)). Since the shear modulus in the crust dominates the magnetic pressure in most parts of the crustal regions for a 
dipole field of strength Bp <10^'G, this treatment may bo justified so long as the crustal toroidal modes arc well decoupled 
from the fluid core. But, if the torsional waves in the crust are strongly coupled with magnetic waves in the core, the treatment 
may not be justified because the angular dependence of the magnetic waves in the fluid core cannot be correctly represented 
by a single spherical harmonics. For example, the analysis by Reese, Fincon, & Rieutord (2004) of toroidal modes in a fluid 
shell have employed flnite series expansions of long length for the perturbations. 

In this paper, using the method of series expansions of perturbations we calculate toroidal modes of a strongly magnetized 
neutron star having a fluid core and a solid crust, where the entire interior is assumed to be threaded by a poloidal magnetic 
field. We employ two different sets of oscillation equations, one for fiuid regions and the other for the solid crust, and solutions 
in the solid and fluid regions are matched at the interfaces between them to obtain an entire solution of a mode. The method 
of calculation we employ is presented in §2, and the numerical results are given in §3, and conclusions are in §4. 



2 METHOD OF SOLUTION 

2.1 Magnetic Fields in the Interior 

We employ an inner magnetic field of Ferraro (1954) type. For a magnetized star in hydrostatic equilibrium, Ferraro (1954) 
assumed an interior poloidal field given by 
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where ?7 is a scalar function, and looked for a particular solution to the partial differential equation: 
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where /t is a constant to be determined. Assuming the scalar function U is given by 

U = (Cir^ + Cir") sm 9, (3) 

and imposing the condition that the interior field is continuous at the stellar surface with the outside dipole field, Ferraro 
(1954) obtained 



Ci = --Bp, 



(4) 



5 _ K _3Bp 

where Bp denotes the strength of the magnetic dipole field at the pole, and R is the radius of the star. For modal analysis 
shown below, we use this equilibrium magnetic field configuration with Bp being a parameter. 



2.2 Oscillation Equations 

We assume the temporal and azimuthal angular dependence of perturbations is given by a sing le factor e*('"'*+"*) for osciUations 
of magnetized, non-rotating stars, where u! is the oscillation frequency in the inertial frame, and m denotes the azimuthal 
wave number. The linearized basic equations governing axisymmetric (m = 0) toroidal modes propagating in the solid crust 

of the magnetized star may be given by 

-a.% = i[V.a']^ + ^[(VxS')xB],, (5) 

S; = [Vx(^xB)]^. (6) 

In equation (5), cr' denotes the Euler perturbation of the stress tensor and is obtained from the Lagrangian perturbation 
defined in Cartesian coordinates by 

6(7ij = {Tipu)Sij +2n{uij - ^u6ij) (7) 
with Uij being the strain tensor defined by 



Axisymmetric toroidal modes of magnetized neutron stars 3 



1 (d^_^d^ 

2 \ dxj dxi 



(8) 



where Sij denotes the Kronecker delta, fj, is the shear modulus, and u = Ym=i ( ^ S ' McDermott et al 1988, and Lee & 
Strohmayor 1996). The linearized equations for toroidal modes in the fluid regions are obtained by simply dropping the term 
p~^V • a' in equation (5). 

Since separation of variables is not possible for oscillations of magnetized stars, we employ series expansions of a finite 
length jmax for the displacement vector ^ and the perturbed magnetic field B' using spherical harmonic functions Yi'"{9,cj)) 
for a given m. We have for axisymmetric toroidal modes 
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where Ij = |m| + 2{j — 1) and 



+ 1 for even modes, and 



\m\ + 2j — 1 and I'j = Ij — 1 for odd modes, respectively. 



and J = 1, 2, 3, • • • , jmax, and Bo is a normalizing constant, which is set equal to Bp. Note that, for axisymmetric toroidal 
modes with m = 0, we have to redefine Ij and I'j as Ij = 2j and I'j — Ij — 1 for even modes, Ij — 2j — 1 and I'j — Ij + 1 for odd 
modes for j = 1, 2, 3, ■ ■ • , jmax. Most of the numerical results shown below are obtained for jmax = 20. In this convention, 
the angular pattern of B'^ (of ^<^) at the stellar surface is symmetric (anti-symmetric) about the equator for even modes and 
it is anti-symmetric (symmetric) for odd modes. 

If we use as dependent variables the vectors t and defined by 



t=(T,,{r)), b^={br.{r)), 



the oscillation equations for fiuid regions are given by 

^Q,Cor-^it = Ciib^ — 
Bo dr 
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and those for the solid crust are given by 

^Q^Cor-^it = Ciib^ - 
Bo dr 
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and G is the gravitational constant, M and R are the gravitational mass and the radius of the star, respectively, and 
Bi=Br/ cos e, B2 = -Be/ sin 6, 



(11) 

(12) 
(13) 

(14) 
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and pn 



Bq/Stt. The definition of the matrices Co, Ci, Qq, Q^, Ao, and Ai is found, for example, in Lee (1993), and the 



matrices Co, Ci, and Ai are defined by 
60 = Co, Ci = {Ci}, Ai=Ai 
for even modes, and 
60 = {Co}, Ci = Ci, Ai = {Ai} 



(19) 



(20) 



for odd modes, where {A} = (Ai,j+i) for a matrix A = {Aij). 

For the inner boundary condition we require that the functions it and ib^ /r are regular at the stellar center. The boundary 
condition at the stellar surface is given by ib^ = 0. The jump conditions at the solid-fluid interfaces is the continuity of the 
function it and of the r4> component of the traction, the latter of whcih leads to 
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where [/(r)]± = lim^^o [fir + e) - /(r - e)]. 



(21) 



2.3 Neutron Star Models 

We integrate the Tolman-Oppenheimer-Volkoff equation (see, e.g., Shapiro & Teukolsky 1983) to obtain zero temperature 
(T = 0) neutron star models, where no effects of the magnetic field on the equilibrium structure are included. The equation 
of state used for the inner crust and the fluid core is that given by Douchin & Haensel (2001). The equation of state in the 
outer crust is that given by Baym, Pethick, & Sutherland (1971), and for the fluid ocean we simply use the equation of state 
for a mixture of a completely degenerate electron gas and a non-degenerate gas of Fe nuclei. The boundary between the fluid 
ocean and the outer crust is set rather arbitrarily at the density p = lO^'g cm~^. 

For the solid crust, we employ the average shear modulus fies (Strohmyer et al 1991), which in the limit of F s 
{Zef/{akBT) — »■ oo is given by 

Mefi = 0.1194 X i^^, (22) 
where n is the number density of the nuclei and a denotes the separation between the nuclei defined by 

= 1, (23) 

and ks is the Boltzmann constant. 

For modal analyses, we calculate neutron star models of mass M = 1.245Mq and M = 0.8967Mq. For the former model, 
the radius and the central pressure and density are R = 1.18 x 10^ cm, Pc = 10'^^ dyne/cm^, and pc = 8.75 x 10^* g/cm'. 
respectively. The ratio of the thickness of the crust, Arcrust, to the radius is Arcrust/-R — 0.093, and the surface ocean above 
the crust is very thin. For the latter model, the parameters are R = 1.19 x 10® cm, = 5 x 10^* dyne/cm^, pc = 6.75 x 10^* 
g/cm^, and Arcrust — 0.139, respectively. 

Figure 1 plots the frequency uj of the crustal toroidal modes of the two models for low values of harmonic degree I, where 
and v = a;/27r, and no effects of magnetic field are included. For each value of i > 2, the lowest frequency 
mode is the fundamental mode that has no radial nodes of the eigenfunction. There exists no fuudameutal mode for I = 1. The 
figures show that the fundamental mode frequency rapidly increases with increasing but the frequencies of the overtones 
remain almost constant for varying I. As first discussed by Hansen k, Cioffi (1980), the frequency of the fundamental toroidal 
mode in the solid crust is rather insensitive to the neutron star mass, and that of the overtones is inversely proportional to 
the thickness Arcrust of the crust, which is consistent with Figure 1. 



3 NUMERICAL RESULTS 

With our numerical method employed to calculate toroidal modes of a magnetized star, we find numerous solutions to the 
oscillation equations for a given Bp. Most of the solutions thus obtained, however, are dependent on jmax, the length of the 
expansions, and we have picked up only the solutions that are independent of jmax. 

In Figure 2 the oscillation frequencies of the toroidal modes of the M = 1.245M0 model are plotted versus the number of 
radial nodes of the expansion coefficient iTij, where the left panel is for Bp = 5 x lO^^G and the right panel for Bp = lO^^G, 
and the filled circles and squares indicate odd and even modes, respectively. The figures clearly show the existence of distinct 
mode sequences, in each of which the frequency of the mode remains rather constant and increases only slowly as the number 
of nodes increases. 

To compare the two cases of Bp = IQ^^G and Bp = 5 x lO^^G, we plot in Figure 3 the frequency ratio aj(10^®G)/aj(5 x 
lO^^G) versus a)(5 x 10^^G)/wo, where the filled squares and circles are for even and odd modes, respectively. Except for the 
lowest frequency sequence of the even modes, the ratio is approximately equal to ~ 2, and it decreases as the radial order 
increases, indicating that the frequency is approximately proportional to the field strength and is largely determined by the 
strength of the magnetic field. This is not the case for the lowest frequency sequence, for which the ratio rapidly increases 
from ~ 1 to ~ 2 as the radial order increases, demonstrating that the response of the oscillation frequency to the field strength 
is not the same as that for the overtone sequences. We note that the lowest frequency mode found in the figure is a much 
more slowly increasing function of Bp than the modes in the overtone sequences. 

It is convenient to write the mode frequency as 



u!p{k, n) = u!b + {k + jp)Aa; + (n — k)5u! + 5{k, n), 



(24) 
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Figure 1. Frequencies of tlie crustal toroidal modes for low degree Is, where luq = GM/R^ and i/ = uj/2tt, and no effects of magnetic 
fields are included to calculate the modes. The left panel is for the M = 0.8967Mq model, and the right panel for the M = 1.245Mq model, 
respectively. For each value of i > 2, the lowest frequency mode is the fundamental mode having no radial nodes of the eigenfunction. 
No fundamental mode exists for I = 1. 
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Figure 2. Frequencies of toroidal modes of the M = 1.245M0 model versus the number of radial nodes of iTi^ , where ojo = 

and 1/ = Ltj/2n, and the filled squares and circles denote even and odd modes, respectively. The left panel is for the case of Bp = 5 x lO^^G, 

and the right panel for Bp = lO^^G. No fundamental modes without radial nodes are found. 



where the non- negative integer k denotes the order of the mode sequences of a given parity (even or odd), and the integer 
n(> k) is the number of nodes of the expansion coefficient iTi, and the integer jp is set equal to for the even mode sequences 
(p = e) and to 1/2 for the odd mode sequences {p = o). Since there exists no fundamental modes that have no radial nodes 
of the expansion coefficients, we require n > 1 even for fc = 0, for which we have to replace n with n — 1 in equation (24). 
The quantities ujb, which is set equal to the lowest frequency obtained, Au), and Su) are assumed independent of k and n, 
and the quantity 6{k, n) is assumed to be substantially small compared with Aw and 6uj. Note also that we find no mode 
sequence having fc = for odd modes, and hence we have to start with fc = 1 for them. Except for the case that involves the 
lowest frequency sequence of even modes, u)p{k + l,n) — uip{k, n) = Aw + 5{k + 1, n) — 5{k, n) does not strongly depend on fc 
and n, and u)p{k, n + 1) — ^^(fc, n) = 5uj + 5{k, n + 1) — 5{k, n) slowly decreases as n increases. It is interesting to note that 
the quantity Aw is approximately proportional to the field strength Sp, but the quantity u)p{k,n + 1) — ijjp{k,n) is rather 
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Figure 3. Frequency ratio a;(10^®G)/a;(5 x lO^^G) versus a;(5 x 10^^G)/a)o for the model of M = 1.245M0, where the filled squares 
and circles axe for even and odd modes, respectively. 
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Figure 4. Same as Figure 2 but for M = O.8967M0. 



insensitive to Bp. For example, for the model of M = 1.245M0, we have Aw/wo — 0.052 for = lO^^G and Acj/tJo — 0.0255 
for Bp = 5 X IQ^'^G, and Slo/ljo ~ 0.0017 for both cases. 

Figure 4 is the same as Figure 2 but for the model of M = O.8967M0 , the crust of which is thicker than that of the 
M = 1.245Mo model. Although the frequency spectra look almost the same between the two models, for a given set of (p, fc, n) 
and Bp, the frequency found for the M = O.8967M0 model is higher than that for the M = 1.245M0 model. This is because 
the fluid core radius R — Arcrust (the Alfven velocity va oc Bp/y^) of the M = O.8967M0 model is smaller (larger) than that 
of the M = 1.245M0 model, and hence the traveling time {R — Arcrust)/«A of magnetic perturbations in the core for the 
former is shorter than for the latter. 

As examples of the cigcnfunctions of the toroidal modes, wo plot the expansion coefficients iT;/ and ihj of three even 

3 

modes as functions of the fractional radius x = r/R'm Figures 5 to 7, where the solid, long-dahsed, short-dashed, and dotted 
lines are the expansion coefficients with j = 1 to 4, respectively. For the lowest frequency sequence even modes plotted in 
Figures 5 and 6, the amplitudes xiTii are much larger than those of ihJ , indicating that the kinetic energy is dominating the 
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0.005515 for Bp = IQi^G, where u)o = \/ GM/R^' and M = 1.245M0. The solid, long-dashed, short-dashed, and dotted lines are for the 
expansion coefBcients associated with j = 1 to 4, respectively. The amplitude normalization is given by iT;/ = 1 at the stellar surface 
X = 1. 
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Figure 6. Same as Figure 5, but for the even mode of aJe(0, 5)/aJo = 0.01671. 



magnetic energy. The magnetic perturbations ibj have amplitudes also in the solid crust. This is not the case for the modes 
in the overtone mode sequences having fc > 1 as exemplified by Figure 7, where the amplitudes of the magnetic perturbations 
arc much larger than xiTii and are strongly confined in the fluid core, having almost no amplitudes in the crust. Also note 
that the expansion coefficients with different values of j seem to coalesce to a single curve in the outer core of the star. 

Figures 8 to 10 give plots of the expansion coefficients zT;/ and ihj for three odd toroidal modes associated with (fe, n) = 
(1, 1), (1,5), and (5,5). The magnetic perturbations ibJ of the modes have much larger amplitudes than xiTii and they are 
strongly confined in the fluid core. The functions xiTii have comparable amplitudes both in the fluid core and in the solid 
crust. As in the case of even modes, the expansion coefficients of the mode (5,5) tend to coalesce to a single curve in the outer 
core region. 
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Figure 8. Expansion coefficients xiT,/ and ibl~ as functions of the fractional radius x = r/R for the odd toroidal modes of uJo(l, ^j/i^o = 

j h 

0.08271 for Bj, = W^G, where wo = ^/gmJb? and M = 1.245M0. The solid, long-dashed, short-dashed, and dotted lines are for the 
expansion coefBcients associated with j = 1 to 4, respectively. The amplitude normalization is given by iT;/ = 1 at the stellar surface 
X = 1. 



We need to examine how the existence of a solid crust affects the frequency spectrum of the toroidal modes. Figure 11 
shows the toroidal mode frequencies of the M = 1.245M0 model for Bp = 5 x lO^^G, where the modes are calculated by 
treating the entire interior as a fluid. Without the crust, the lowest frequency sequence is composed of odd modes, and the 
even mode sequence corresponding to A; = 0, which appears as the lowest sequence when the solid crust is included, is missing. 
Therefore, for the fluid star we have to set jp equal to jp = for odd mode sequences with p = o and to jp = 1/2 for even 
mode sequences with p = e in equation (24) . In the lowest frequency odd mode sequence, there exists the fundamental mode 
having no radial nodes of the expansion coefficients, corresponding to {k,n) = (0,0). Since the oscillation equations (12) and 
(13) for fluid regions contain the oscillation frequency only in the form of the ratio uj/Bp, the frequency of a given (p, k, n) is 
exactly proportional to the field strength Bp. This suggests that the quantities ujb, Au>, Suj, and 6{k,n) defined in equation 
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Figure 10. Same as Figure 8, but for the odd mode of Wo{5, 5)/wo = 0.2976. 

(24) are also proportional to Bp. Note that if we write ojf, oc for the lowest frequency mode for the case with a solid 
crust, where the exponent /? > may be weakly dependent on Bp, we have (3 <l/4 for Bp >10^^G, the feature of which is 
different from the case without the crust. Figure 11 indicates that the frequency separation between the lowest and second 
lowest sequences is approximately ~ Aw/2 for the case without the crust, but this separation is given by ~ Auj for the case 
with the crust. Note that the frequency separation between the second and third lowest sequences, for example, is given by 
~ Aaj/2 for both cases. 

As suggested by Figures 2, 4, and 11, it is not always possible to numerically find jmax independent modes associated 
with large values of n and k. Possible numerical reasons for the difficulty may be that the modes we want axe immersed in a 
dense spectrum of j'max dependent solutions, particularly for weak Bp, and that a large jmax and high numerical precision are 
required to correctly calculate the expansion coefficients in the vicinity of the stellar center for modes associated with large 
n> k. Prom the theoretical point of view, however, it is not necessarily clear that the toroidal modes of the kind we find in 
this paper exist for arbitrary sets of integers A;(> 0) and n(> k). 
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respectively. The modes are calculated by treating the entire interior as a fluid. 



4 CONCLUSIONS 

In this paper, we have calculated toroidal modes of magnetized stars with a solid crust, where the entire interior of the star 
is assumed to be threaded by a poloidal magnetic field that is continuous at the stellar surface to the outside dipole field. We 
find distinct mode sequences of the toroidal modes, in each of which the mode frequency remains rather constant and only 
slowly increases as the radial order of the modes increases. In the presence of a solid crust, the frequency separation between 
the lowest and second lowest frequency mode sequences is approximately given by Aw, but that between the second and third 
lowest frequency mode sequences by Acij/2, where the frequency separation Ao; is roughly proportional to the field strength 
Bp. This frequency pattern of the low frequency sequences is different from that found for the model without the crust, for 
which the frequency separation between the sequential sequences of low frequency modes is given by Aa'/2. We also find that 
for the equation of state we use, /S.lo is larger for smaller M. 

The eigenfunction of the modes belonging to the lowest frequency sequence have much larger amplitudes than B'^ and 
can penetrate into the solid crust. On the other hand, (,4, of the modes belonging to the higher frequency sequences is much 
smaller than B'^, which is well confined into the fiuid core and does not have any substantial amplitudes in the solid crust. 

The frequency ranges of the toroidal modes we find for the magnetized neutron star with Bp ~ 10^'' G overlap the QPO 
frequencies found for the magnetar candidates, SGR 1806-20 and SGR 1900+14. This suggests that we may interpret the 
observed QPOs based on the magnetic toroidal modes, and that detailed comparisons between observed frequency spectra 
and theoretical calculations make it relevant to infer physical parameters of the magnetar candidates, such as the equation 
of state and the strength of the magnetic field. But, we think it worth pointing out that except for the modes belonging to 
the lowest frequency sequence, the magnetic perturbations, which have much larger amplitude than are well confined in 
the fluid core and do not have substantial amplitudes in the solid crust, which make it difficult for the modes to be directly 
observable. If the magnetar candidates do no have a crust, the problem of observability could be avoided. In this case, however, 
we have to use more realistic surface boundary conditions than B' = used in the present calculations. Note that, although 
the existence of a solid crust affects the frequency pattern of the low frequency mode sequences, the frequency range of the 
magnetic modes itself is not very much dependent on the presence or absence of a solid crust unless the crust is extremely 
thick. 

We have tried to find toroidal modes well confined in the solid crust or core toroidal modes that are in resonance with the 
crustal toroidal modes, but failed. We also find it extremely difficult to identify distinct toroidal mode sequences for magnetic 
fields weaker than Bp <10^^G when a solid crust is included in the models. It is therefore not clear whether the toroidal 
mode sequences of the kind we find in the present paper can survive also for weakly magnetized neutron stars with a solid 
crust. If the field strength is much weaJier than Bp ~ lO^^G, the magnetic fields may have only minor effects on the crustal 
toroidal modes (e.g., Lee 2007), and it will be justified to use frequency spectra of the crust modes theoretically obtained 
in the weak field limit to interpret observed QPOs. As briefiy noted in the last section, it is not theoretically clear how the 
frequency spectra of the toroidal modes of magnetized stars should look like, and how the toroidal modes behave in the limit 
of n > A; — > 00. We may even speculate that the frequency spectra we obtained reflect the existence of continuous spectra (see, 
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e.g., Goedbloed & Poedts 2004, see also Levin 2007), but the detailed numerical analysis of continuous frequency spectra is 
beyond the scope of this paper. It will be worthwhile to examine the effects of an interior toroidal field on magnetic modes. 
It is also needed to extend the present analysis to a general relativistic formulation (e.g., Sotani et al 2006, 2007). 
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